function [Z ] = isStableSys( dx,dy,u )
%ISSTABLESYS Summary of this function goes here
%   Detailed explanation goes here
param.a=1.446;
param.b=1.408;
param.c=1.4370;
param.I=1549.034;
param.m=2220;
param.g=9.81;

for i=1:size(dx,1)
    for j=1:size(dx,2)
        x=[dy(i,j);dx(i,j);zeros(4,1)];
        [A,~,~,~]=carLin(x,u,param);
        Z(i,j)=double(~any(eig(A)>0));
    end
end
end

